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We show that extreme vacuum pressures can be measured with current technology by detecting 
the photons produced by the relativistic Thomson scattering of ultra-intense laser light by the 
electrons of the medium. We compute the amount of radiation scattered at different frequencies 
and angles and design strategies for the efficient measurement of pressure. In particular, we show 
that a single day experiment at a high repetition rate Petawatt laser facility such as VEGA, that will 
be operating in 2014 in Salamanca, will be sensitive, in principle, to pressures p as low as 10" 16 Pa, 
and will be able to provide highly reliable measurements for p > 10 _14 Pa. 

PACS numbers: 42.62.-b, 07.30. Dz, 41.60.-m, 52.38.-r 
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Introduction.- Pressures corresponding to Extreme- 
High Vacuum (XHV), p < 10 _10 Pa [1], are measured 
by ionization methods: the atoms in the sample are ion- 
ized and the produced charged particles are collected by 
applying an electric field. This procedure is fully reliable 
for pressures as low as 10 _11 Pa @|. Although there are 
techniques able to push this limit down|3[, its use would 
be questionable since in this regime the electron stimu- 
lated desorption, the so-called X-ray limit, or the out- 
gassing from the hot cathode cannot be neglected 0, Q • 
It is therefore of crucial importance the introduction of 
a new, alternative method aimed at providing an inde- 
pendent measurement of XHV pressures below 10 _11 Pa 
without significantly altering the pressure itself and free 
of the aforementioned limitations. 

In this Letter, we propose the idea of using photons to 
gauge the extreme vacuum properties. For this purpose, 
the advent of ultra-high intensity lasers [H[ has provided 
a new class of light sources which are powerful enough to 
produce a measurable signal even in conditions of XHV. 
In the XHV, the remnant pressure is essentially produced 
by the hydrogen released by the walls. When interact- 
ing with high-intensity laser light, the electrons can be 
considered as free, therefore the main source of dispersed 
light is relativistic nonlinear Thomson scattering. This 
process was studied in detail in (il-Hoj and experimentally 
observed in [ll|, and ma y b e used to measure the peak 
intensity of a laser pulse [12j . 

We will compute the number of scattered photons as 
a function of the electron density of the medium and the 
parameters of the ultra-intense laser pulse (wavelength, 
peak power, waist radius, pulse duration and repetition 
rate). In the XHV regime, collective effects of the elec- 
trons as those discussed in [H, ^ can be neglected. The 
number of scattered photons is proportional to the num- 
ber of scattering centers, which is proportional to the 
pressure. Harmonics are generated in the scattering pro- 
cess. It will be shown that most of the scattered photons 
correspond to the incident wavelength (n = 1). Never- 
theless, detection of photons with n = 2,3, .. . may be 
possible and useful. Remarkably, we find that it should 



be possible to provide a highly reliable measurement of 
pressures as low as p « 10 -14 Pa in realistic conditions 
at facilities that will be available in the near future. As a 
consequence, ultra-intense lasers may be able to push the 
physical limits for measuring a basic magnitude like pres- 
sure. A lower cost practical application of this result can 
be the use of more common intermediate-intensity lasers 
as an alternative instrument of measuring high vacuum 
pressure in a non-extreme regime. 

On the other hand, ultra-intense lasers need high vac- 
uum to operate and it is important to determine their 
operation conditions. Moreover, extreme vacuum is a 
necessary requirement of many of the experiments based 
in ultra-intense lasers that have been proposed in the last 
few years aimed at demonstrating the quantum vacuum 
polarization [l3l - [l5| and at searching for new particles (Tfj. 
It is fascinating in these cases that the laser itself may 
provide an efficient tool to monitor the pressure in the 
chamber, substituting or complementing other conven- 
tional methods. 

Relativistic Thomson scattering.- Our computations 
are based on some of the results of |7| , which we briefly re- 
view for completeness and to fix notation. We introduce 
the dimensionless parameter related to the intensity / 
and wavelength Aq of the beam as: 
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where ro ~ 2.82 x 10 _15 m is the classical electron radius. 
Let us also define A4 = 1 + -|g 2 sin 2 (#/2). Relativis- 
tic effects play a role for q > 1, corresponding to I > 
2 x 10 18 W/cm 2 (for Ao = 800nm). When a linearly po- 
larized plane wave impinges on a free electron, the power 



scattered per unit solid angle is: 
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n is the harmonic number and is dimensionless: 
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with cos a = sin 6 cos(/?, where 6 G [0,7r], (/? G [0, 27r) 
are usual spherical coordinates. Forward scattering cor- 
responds to = and ip = 0, 7r point along the polar- 
ization axis. The F s n can be written in terms of Bessel 

functions as: F s » = £+ = -oo * ( ^Tm^ ) * 
[J 2l+n+s (^P) + J 2i+ „- s (^P)]. These results 
hold in the laboratory frame, in which the scattered wave- 
length is shifted as A^ n ^ = MXo/n. 

In this paper, we will only consider linearly polarized 
laser pulses. The main difference in the case of circular 
polarization would be the lack of the dependence on the 
azimuthal angle. 

Modeling a realistic situation.- Our main goal is to 
compute the average number of photons scattered when 
a laser pulse traverses a vacuum chamber. It is cru- 
cial to take into account the distribution in space of 
the incident radiation. We model the pulse as having 
a Gaussian profile (waist radius wq). The intensity then 

/ \ 2 2r 2 

reads: I = Io ( ) e w(z)2 • The beam radius evolves 

as w(z) = woy/l-\- z 2 1 z\ where the Rayleigh range is 
zr = ttwq/X. Our estimate for the scattered radiation 
will be obtained by performing the appropriate integral 
after inserting the intensity profile in (j2j). In particular, 
the number of photons of the n'th harmonic produced by 
a single pulse are = n e r J ^^^^dQd 3 x, where 

n e is the number of electrons per unit volume, which we 
will assume to be uniform and r the pulse duration. It is 
useful to rewrite the integral in terms of the dimension- 
less quantities p = r /wo, £ = zj zr, such that: 

where go is related to the intensity at the beam focus. 
After some simple manipulations one can write: 

N^=K J J P rW(q)dpdt, (4) 

where we have introduced a function of q: 

Y^ n \q)= / -f^MsmOdOdip (5) 
Jo Jo n 

and the parameter: 

1 ,tt 2 ^ 
JC = -n e (cr)— -2-a, (6) 
z a 

where a « is the fine structure constant. Notice 
that the integral in Eq. ^ only depends on n and on 
qo whereas the rest of quantities describing the physical 
situation are factored out in JC. 

We have computed the number, frequency and spatial 
distribution of the photons that may be produced as a 
function of the incoming laser pulse parameters, by nu- 
merical integration of the expressions in (J2J), (|5]). Fig. 
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[H shows a plot of the function T n (q) for n = 1,...,4. 
Since it is impossible to have a detector covering the 
full solid angle, we also show the results when the in- 
tegral is performed over a reduced range of the polar 
angle G (0 cut ,7v — cut ). In any realistic situation one 
should have 6 cut ^> 0d, being 6 a = Xo/(ttwo) the beam 
divergence. 
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FIG. 1: (Color online) In solid lines, the function r (n) (q) 
found by numerical integration. Dashed lines are found by 
cutting the ^-integration in Eq. © with cut =0.1. 

Using these results for Y^ n \q) and the Gaussian inten- 
sity distribution ((3]) one can readily compute the integral 
in (j4]). We have confined the integration to the region 
where an electron of a hydrogen atom can be considered 
as free q > q cut = 0.01, see discussions below for more 
details. The results are plotted in Fig. [2l We find the 
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FIG. 2: (Color online) Semi-logarithmic plot of N^ n) /(Kql) 
for n = 1, . . . , 4. 

following asymptotic behavior (valid for qo large enough 
depending on n): 

N^^c n Kql (J) 

with the values a w 275, c 2 « 1.3, c 3 w 0.22, c 4 « 0.088. 
If for Fig. [2] one performs the same cut as before in the 
integration region 0.1 < < tt — 0.1, the correction to the 
result is tiny, below 1%. This happens because most of 
the photons are not generated at the maximum intensity 
region — the beam focus — , but at the larger volume 
where the Gaussian profile presents moderate values of 
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q, irrespective of how large qo might be. In particular, 
for n = 1 most of the photons are generated in a re- 
gion with q < 1 and therefore one can find an approx- 
imation to the n = 1 result using the simpler expres- 
sions for non-relativistic Thomson scattering. One gets 
ci w 87r/(9q cut ). 

In order to understand qualitatively the results for 
n > 1, we may approximate the plateaus of T^^q) dis- 
played in figure [1] by Heaviside step functions T^ n \q) = 
b n Q(q — q s tep,n)- Then, defining the limits of the q > 
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1, one can estimate the integral in Eq. 
p b n dpd^ ~ q b ^ q ° where we have only 
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kept the leading term in qo/q s tep,n- This simplified anal- 
ysis fits qualitatively the numerical results and explains 
the cubic dependence in qo: it is a consequence of the 
fact that, roughly, both pu m , £u m grow linearly with go- 
In Fig. [3j we present plots of the value of r( n ) (q(p, £)), 
which show the region in space in which the photons are 
scattered. The information displayed in Fig. [3] indicates 
which the precise region of the chamber where the elec- 
tron density measurement is taking place. 
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FIG. 3: (Color online) T (n) (<?(>, £)) in two cases. For n > 1, 
photons are produced in the region where q > 1. The region 
where n — 1 photons are scattered is larger since they are the 
only outcome of the non-relativistic regime q <C 1. 

In Fig. |4]we plot two examples of the angular distribu- 
tion of the emitted photons, obtained by performing the 
integral in dU on the p — £ space and leaving it as a func- 
tion of and (p. For a given harmonic, the angular depen- 
dence does not change too much when modifying qo. The 
reason for this is that — as noted above — even when 
qo is large, a copious amount of radiation comes from 
the region of smaller q. This same argument explains 
why the distributions are not forward peaked, as one may 
naively expect. In fact, the plot for n = 1 can be hardly 
distinguished from the angular distribution correspond- 
ing to linear Thomson scattering sin 0(1 — sin 2 #cos 2 <p). 
These considerations may be useful when looking for an 
optimised configuration of photon detectors in an actual 
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FIG. 4: (Color online) Plots of 1/n f f pf {n) M sm6dpd£ as 
a function of and <p, see Eqs. Q, J5}. 

We now turn to the frequency distribution of 
the scattered photons. Harmonics are not emitted 
in multiples of the original frequency but there is 
a shift related to electron recoil A^ n ^ = A4\o/n, 
where M. depends on q and 0. Formally, we can 

find the spectral distribution by writing d} J n) = 
KjdpjdijdOjdif (p±fMMsm06(\M - ^)). 
In Fig. [5j we plot the result of the numeric integration 
for qo = 3. The result does not appreciably change for 
larger values of qo- In the final part of this Letter, we 
will comment on limitations of this computation. 
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FIG. 5: (Color online) Differential frequency distribution of 
the scattered radiation for n — 1, . . . , 4, in logarithmic scale. 



Quantitative estimates.- Since the number of scattered 
photons is directly proportional to the number of scat- 
tering centers, at extreme low pressure, the scattered ra- 
diation will be extremely weak. In order to get an idea 
of whether a given pressure can be measured with this 
method, we must know the number of photons that can 
be detected in a reasonable amount of time. If we assume 
that only atomic hydrogen is left in the chamber, the re- 
lation between the electron density and the pressure is 
n e = p/(kBT). Using the result (pp), we can estimate the 
total number of photons detected in a period of time At 
in terms of the laser repetition rate r r , the total energy of 
each pulse E pu i se = r Iqtt Wq/^ and the detector efficiency 
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/ — which include geometric and quantum factors. 

It might seem weird that for fixed energy the signal grows 
with the waist radius. What happens is that having 
smaller maximum intensity is compensated by a larger 
interaction region. Nevertheless, when wq it too large, 
qo becomes small and ([7j) and © lose their validity (see 
Fig. EJ. 

Eq. dSJ is the main result of this paper. To be con- 
crete, we can now evaluate the limiting pressure that can 
be measured with this method at a given ultra-intense 
laser facility in a reasonable time span. As an example, 
we will consider the Petawatt laser VEGA that will be 
available in 2014 at the CLPU of Salamanca [ItJ, hav- 
ing repetition rate as large as r r = Is -1 , with pulses of 
Ao = 800nm, r = 30fs and E pu i se = 30J. Taking e.g. 
T = 300K, wo = 20/im and a day run, At = 1 day, 
and assuming an efficiency / = 0.5 for n = 1 which is 
a realistic value for commercially available single photon 
detectors at Ai = 800nm, we can compute the limiting 
pressure that can be measured within 3 standard devi- 
ations by taking N^ det = 10 in Eq. (jSJ). We obtain 
Plimit — 10 _16 Pa. Of course, this sensitivity should be 
corrected by a geometric efficiency factor, depending on 
the effective area of the detector that is chosen. Geomet- 
ric efficiency corrections will be more important for the 
n = 1 photons since the region where they are scattered 
is large, Fig. 02 Note that the angular cut that we have 
imposed in our computation ensures that the detected 
photons will not be confused with those of the beam that 
do not undergo Thomson scattering, that would give no 
observable signal in the integration area for one day run 
(other sources of noise will be considered below). An- 
other way of avoiding such kind of background would be 
the measurement of the n = 2 harmonic. From Eq. flSJ, 
we obtain that the limiting pressure that can be mea- 

(2) 

sured by detecting det = 10 photons after one day 
run would be pnmit — 10 _14 Pa, if all the other parame- 
ters are taken as above except the efficiency, that can be 
as large as / ~ 0.6 for A2 — 400nm in state-of-art sin- 
gle photon detectors. For all these reasons, we conclude 
that the detection of the n = 2 harmonic will provide 
an independent measurement of the pressure above the 
10 _14 Pa range for a 1 day run at VEGA, complemen- 
tary to the more sensitive measurement due to the n = 1 
wave. Taken together, these measurements could be used 
for a kind of self-calibration of the whole procedure. The 
result for the measurement of pressure as low as 10 ~ 14 Pa 
would then be highly reliable, provided that the detec- 
tion is accurate enough and the noise level can be kept 
below the signal, which are feasible tasks with present 
technology as we discuss below. 

Background analysis.- A source of noise are thermal 



photons, whose expected number per shot is 

N*% t ~ 2tt cS det At { dX f{X) f 4 (9) 

where /(A) is the sensitivity of the detector as a func- 
tion of the wavelength A, Sdet is the detecting area and 
T is the temperature of the vacuum tube. At ordinary 
temperatures T ~ 300if , this background can be made 
completely negligible in all the configurations that are of 
interest for the present work by using a wavelength fil- 
ter on the detector, cutting off all the wavelengths larger 
than Ao + AAo, where AAo — Xq/(ct) is the uncertainty 
in the pulse wavelength. 

A potentially higher source of noise is due to the dark 
counts of the detector, that can be kept below 10s _1 in 
avalanche photodiodes featuring high efficiencies at the 
wavelengths discussed in this Letter. If we require that 
during each repetition the detection window is opened for 
a very short time, which can be as short as two nanosec- 
onds with present technology, we can ensure that after 
the ~ 10 5 repetitions in the 1 day experiment at VEGA 
the total dark count would be unobservable. Such gating 
of the detector would also provide an efficient protection 
mechanism against the backscattered photons from the 
walls of the vacuum tube. In fact, assuming space dimen- 
sions of the tube of the order of few tens of centimeters or 
larger, the backscattered photons would reach the detec- 
tor out of the detection window. A promising alternative 
could also be the use of superconducting single photon 
detectors [HI [l9[, that are able to reduce the dark counts 
below 10 _2 5 _1 in both the visible and infrared ranges. 

Validity of approximations.- We discuss now several 
approximations and assumptions that have been made 
in deriving our results. First of all, radiation reaction 
and quantum effects have been neglected in Eq. ((2]), 
which is a good approximation since q 2 <C Ao/ro and 
nhc/Xo <C m e c 2 in relevant situations. Moreover, ((2]) 
is valid for a plane wave. This means that the Gaus- 
sian beam radius should be larger than the transverse 
displacement of the electron which is typically of or- 
der q A. Namely, the formalism is valid for wq ^> qo^o 
and cannot be used for diffraction limited beams. We 
have considered the radiating electrons as free. This is 
valid as long as the atomic potentials can be neglected 
in the presence of the laser beam, namely in the bar- 
rier suppression regime [20|, which for hydrogen corre- 
sponds to I > I B s = 1.4 x 10 14 W/cm 2 [2l|. Since this 
limiting value corresponds to q <C 1 (it is q « 0.01 for 
Ao = 800nm) — and q ~ 1 is related to the onset of rela- 
tivistic effects — , we conclude that the binding energy of 
the electrons does not play a role in the harmonic gener- 
ation we have discussed. For the same reason, harmonic 
generation coming from electron-proton recombination is 
suppressed (and would be further suppressed if polariza- 
tion is non-linear) leaving relativistic Thomson scattering 
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as the dominant process. It should also be mentioned 
that we have used expressions for electrons initially at 
rest since at room temperature they are far from being 
relativistic. 

Finally, we have made a rough modeling considering 
square pulses and not taking into account the effects of 
the time envelope of the pulses. For short pulses — with a 
small or moderate number of light cycles — this approx- 
imation is not accurate, as it has been discussed both 
and quantum [24] frameworks. The 
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main correction that will appear is a spectral broadening 
- called ponderomotive broadening in [22]. Because of 
this phenomenon and of the fact that the initial beam is 
not monochromatic, the plots of Fig. [5] might markedly 
underestimate the spectral width of the produced har- 
monics for short pulses. Corrections may also multiply 
Eq. (j8j) by a factor of order 1, but we do not expect them 
to change the orders of magnitude. However, it could be 
pertinent to make a full study when dealing with a par- 
ticular situation. 

Conclusions.- We have computed the amount and 
spectral distribution of the photons that are produced 
when a gaussian laser pulse crosses a vacuum tube. With 
present detector and ultra- intense laser technologies, this 
implies the possibility of measuring pressures as small 
as 10 _16 Pa. This technique can be self-calibrated and 
highly reliable above the 10 _14 Pa scale. 
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